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Abstract. In this paper we prove, extend and review possible mappings between 
the two-dimensional Cluster state. Wen's model, the two-dimensional Ising chain and 
Kitaev's toric code model. We introduce a two-dimensional duality transformation 
to map the two-dimensional lattice cluster state into the topologically-ordered Wen 
model. Then, we subsequently investigates how this mapping could be achieved 
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physically, which allows us to discuss the rate at which a topologically ordered system 
can be achieved. Next, using a lattice fcrmionization method, Wen's model is mapped 
into a series of one-dimensional Ising interactions. Considering the boundary terms 
with this mapping then reveals how the Ising chains interact with one another. 
The relationships discussed in this paper allow us to consider these models from 
two different perspectives: From the perspective of condensed matter physics these 
mappings allow us to learn more about the relation between the ground state properties 
of the four different models, such as their entanglement or topological structure. On 
the other hand, we take the duality of these models as a starting point to address 
questions related to the universality of their ground states for quantum computation. 
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1. Introduction 

Studying tlie pliysical properties of lattice spin systems is one of the most challenging, 
yet most fascinating areas of condensed matter physics. These systems exhibit a plethora 
of exotic quantum phases, some of which are even believed to be useful for quantum 
computation. Those proposals have in common that they exploit special ground state 
properties of a spin system, e.g. its entanglement or topological structure, to overcome 
the notorious problems of the standard circuit model of quantum computation, such 
as the individual control of qubits or decoherence. Two prominent examples of such 
proposals are Measurement-based Quantum Computation (MBQC) and Topological 
Quantum Computation (TQC). 

In the case of MBQC, the computational task is achieved via single qubit 
measurements on a highly-entangled state [1]. These measurements are not only much 
simpler to realize experimentally, but this approach also reduces the problem of quantum 
control, because it is not necessary to interact qubits after the state is prepared. If the 
entangled state is appropriately chosen, MBQC is equivalent to the standard circuit 
model, i.e. it provides us with a universal quantum computer. One example of such a 
state in two dimensions is the so-called Cluster State [2]. 

TQC, on the other hand, uses a two-dimensional topologically ordered system with 
non-Abelian braiding statistics, and braids the quasi-particle excitations [3] in such a 
way that logical qubits are processed to perform computational tasks [1]. This model 
of quantum computation is robust against errors because the long-range entanglement 
of the system prevents perturbations of small areas of physical qubits from causing 
errors to the logical qubit. A simple model with topological order and Abelian braiding 
statistics which may help us better understand TQC is Kitaev's toric code model [SIE], 
which was originally designed to obtain a topologically protected quantum memory. 

While these two proposals make use of two very different properties of the systems 
to achieve quantum computation, they share the rare property that they are the ground 
states of some non-trivially interacting exactly-solvable spin system. The present work 
shows a method of directly mapping between the two-dimensional cluster state on a 
square lattice and Wen's model, which in turn relates the cluster state with the toric 
code [7j and the Ising model [Hj. In that case, the boundary terms transform nontrivially 
and the ground state degeneracy is changed. These mappings are summarized in 
figure [H This work can be considered a generalisation of the one-dimensional work 
done in reference [9], which transforms the one-dimensional Cluster State into the 
one-dimensional Ising model. This particular work considers the boundary terms and 
topological properties of these models. 

Motivated by the importance of two-dimensional spin models for the realization of 
a quantum computing device, and fascinated by the possibility to map these models 
which seem unrelated at a first glance onto each other, the goal of this work is to review 
and extend known results as well as to prove new possible mappings between these 
systems using our new mapping between the cluster state and Wen's model, and some 
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Figure 1. Summary of the mapping between the two-dimensional cluster state, the 
Wen model, the Ising model and Kitaev's toric code model. The two-dimensional 
cluster state is mapped up to the boundary terms to Wen's model using a duality 
transformation. Wen's model can be diagonialized using a fermionization approach 
leading, up to boundary terms, to a Hamiltonian of anti-ferromagnetic Ising spins. On 
the other hand, the Wen model can be mapped exactly into Kitaev's toric code model 
when the lattice is of even by even dimensions. 



existing mappings. 

This paper is organized as follows: We start with a review of the cluster state model 
and study its topological properties. Next, we prove that in the thermodynamic limit 
the two-dimensional cluster state model can be mapped into the topologically ordered 
Wen's model [10] using a global unitary transformation. In the case of open boundary 
conditions, Chen and Hu [8] have shown how Wen's model can be diagonalized using a 
Jordan- Wigner transformation pTl [T2] . We first review their technique and then turn 
our attention to the more involved case of periodic boundary conditions. We do not only 
give a constructive way how to obtain the spectrum of the Wen Hamiltonian in this way, 
but we also include a discussion about the ground state degeneracy. Finally, we review 
the local unitary transformation that takes Wen's model into Kitaev's toric code model 
[7], establishing a link between the two-dimensional cluster state and Kitaev's toric code. 
This relationship is then used to interpret the boundary terms of the mapping between 
the cluster state and Wen's model when the cluster state on a square periodic lattice is 
transformed. 

2. The Family of Cluster States 

Consider a D-dimensional square lattice and associate to each site a spin-half particle. 
Then, the cluster state, also known as a graph state when it's on a square lattice, on 
that lattice is defined as the common eigenstate of the "stabilizers" Si = Xi Y{j£j\f(i) Zj 
with eigenvalue -|-1 Vz. Here, Xi and Zi are the Pauli spin operators acting on the i-th 
site and M{i) denotes the neighboring sites of the i-th spin. Any cluster state is the 
non-degenerate ground state of the Hamiltonian He = —J2iSi with eigenvalue — A^. 

More insight into the structure of the family of cluster states can be obtained using 
its operational construction: A cluster state can prepared via an application of controlled 
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phase gates Uij = [1 + Zi + Zj — ZiZj]/2 to the product state |+)®^ in the following 
way: 

\c) = {l[u^,m) i+r"" ■ (1) 

The operational representation allows us to write down immediately the complete 
spectrum of He'. If we define \Ch) = Y[j{Zj)"'^\C), where n = {ni,n2, ■ ■ ■ yU^) and nj G 
{0, 1}, then one can easily show that Si\Cfi) = (— 1)"'|C^) Wi and that (C^|C^) = 6fi^rn 
holds. Further, we can read off that the energy eigenvalue of \Cft) is a function of the 
number of spin flips, Eft = —N + 2Y,ini. 

The cluster state is a highly entangled state in a sense that the state has the largest 
relative entropy that any deterministic state can reach for a given number of qubits 
[T3] . The unit of cluster state entanglement is N/2 in terms of the relative entropy of 
entanglement [Hj and it scales like in other measures of multipartite entanglement 
as well [T5] . 

Next, we discuss the topological properties of the family of cluster states defined 
in Eq. ([T]). As it has been proposed in reference [16] and [IT], the von Neuman entropy 
serves as a measure to characterize the topological order of a many-body system 5* in 
the following way: Let A be a subsystem of S, and let L = \dA\ be the length of the 
boundary between A and S* \ A. Then, the von- Neuman entropy is given by 

SM) = -Mpa log pa) = - 7 + ■ ■ -, (2) 

where pA denotes the reduced density matrix of the subsystem A. The ellipsis represent 
terms that vanish in the limit of L — )■ cxd. The scale invariant part 7 characterizes the 
global feature of the entanglement in the ground state, called topological entanglement 
entropy. 

Making use of the operational representation of the cluster state given in Eq. ([1]), 
the topological entanglement entropy of the cluster state can be readily calculated. To 
this end, we make use of the fact that any cluster state can be obtained via a local 
application of controlled phase-gates U^^ and that every individual spin of the system 
is in a completely mixed state. If we further denote as A/" those neighboring sites of the 
boundary OA which lie outside of region A we have that pA = TrA^ UU^ ^\C)a{C\W\. 
Here, \C)a denotes the cluster state defined on subsystem A alone, lj\f is the state of the 
spins in boundary M and U = (JlkedA U'k) is the unitary of the controlled phase gates 
applied to the spins on dA and A/". Thus, pi = ■^Y.{nkedA}=o ^n\C) i{C\Bn, where 
Bfi = UkedAi^k)"''' , n = {ni,n2, ■ ■ ■), and L = Y^kedA is the number of qubits in the 
boundary dA. In this case, the boundary dA is defined. With this result at hand, we 
use Eq. (|2]) to arrive at 

5(|C)) = -^log^ E =L. (3) 

{n^G boundary}=0 

A comparison with Eq. ([2]) reveals that the topological entanglement entropy 7 = 0. 
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a. b. 




Figure 2. a. Shows the standard coordinates, the arrowheads show the direction in 
which the coordinates increase, b. Shows the new coordinates on a 6 x 6 lattice. The 
red arrows shows the j' coordinates still increase, but the diagonal blue lines don't have 
arrowheads because they show lines under which the i' coordinates remain constant. 

Thus, a cluster state on an arbitrary graph is the unique ground state of a stabihzer 
Hamiltonian without topological order. However, as we will show in the next section, 
the two-dimensional cluster state is linked to topologically ordered systems via a simple 
duality transformation. This seems like a contradiction because the energy spectrum 
will not be conserved between a cluster-state model with a non-degenerate ground state 
and a topologically ordered model. This is resolved later in the paper where we consider 
boundary terms. 

3. Mapping of the Two-Dimensional Cluster State to a Topologically 
Ordered System 

In this section, we show how the two-dimensional cluster state Hamiltonian can 
be transformed into a Hamiltonian with topological order using a simple unitary 
transformation. Recall that the Hamiltonian of the two-dimensional cluster state on 
a square lattice is of the form 

The duality mapping between the cluster-state model and Wen's model can be written 
as follows: Firstly, the indices are transformed {i,j) {i'lj') such that j = j' and 
i' = i — j + 1. The effect of this index transformation is shown in figure 2. Then the 
mapping can be conveniently written 

and 

k=l 
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Figure 3. Graphical representation of one term ^Z^-i jZi+i.jZ; j+i of the 
2D Cluster Hamiltonian. 



where the modes are ordered in a diagonal way on the lattice as depicted in Fig |3l 

A straightforward calculation shows that under this dual transformation the cluster 
state Hamiltonian is transformed into Wen's model [TO] 

= - 5Z f^lj-lf^i-l,jf^i-l,j-lf^lj^ (5) 

where we have now reverted back to the original indices. 

Wen's model is of special interest since it exhibits a new kind of phase transition 
where a change of the sign in the Hamiltonian leads to a change of the quantum order of 
the system while the symmetry is preserved. The model is different from Kitaev's toric 
code model [5], but it can be mapped into it through a local unitary transformation. 
Further, Wen's model corresponds to the low energy sector of Kitaev's honeycomb lattice 
model in the limit Jx, Jy P^ - 

This mapping can be discussed in more depth by considering the physical operation 
that performs such a transformation. The mapping is done using the following circuit 
[T8] along each of the diagonals of a periodic lattice oi N x N spins. 

N-l 1'^ 

n {h',N^j ® |0) {0\,^^_j^, + Xe,N-j ® |1) (l|i',iv^,+i; 



N 



The transpose acting on the product of the controlled-not gates reverses their order (of 
course the controlled-not gate is symmetric, so all the transpose does is reverse the order 
of the individual gates). The order is important because all the controlled-not gates have 
to be performed sequentially. They cannot be performed simultaneously because they 
don't commute with one another on pairs of nearest neighbour sites. With this in mind, 
and assuming that each of the controlled-not gates requires only one unit of time to be 
performed, and all the Hadamard gates can be performed simultaneously (as they all 
occur on different sites) and each of the t/j/ operators can be performed simultaneously, 
then the time it takes to transform a cluster state into a topologically ordered system 
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scales linearly with the length of the boundary of the lattice. This agrees with the 
result found in reference [19] who find that the optimal time to generate a topologically 
ordered system scales with the boundary of the system, and not with the volume of the 
system. This matches the rate that topological order is generated in reference [6] , which 
is also mentioned in reference |19j . 

4. Wen's Model and its explicit solution 

Wen's model describes spins on a two-dimensional lattice which are subject to spin-spin 
interactions on 2 x 2-plaquettes [TO] . 

H^ = gY,W,„ (6) 

id 

where Wij = XjjZj+ijZjj+iXj+ij+i and g is the interaction strength. This model 
was designed to be exactly solvable using a projective construction, which is the same 
method adopted by Kitaev to solve the honeycomb model [17]. This paper uses an 
alternative approach to diagonalize this Hamiltonian introduced by Chen and Hu in 
Ref [8] and generalizes this approach to periodic boundary conditions. This method 
uses the Jordan- Wigner transformation to write the spin-half operators in terms of 
fermionic creation and annihilation operators 

H- = ( n n^^Y ) iiiY^'A (^..+^%), (7) 
\j'<j i' J \i'<i J 

that obey canonical anti commutation relations (CAR) {c[.,q} = 6ki. Equivalently, one 
can use the hermitian Majorana operators 

that also fulfill the CAR. Now we define new fermionic modes as linear superpositions 
of Majorana fermions on two neighboring sites {i,j) — {i,j + 1), dij = (ajj + i/3jj+i)/2. 
With these new modes, it follows that 

ioiij/3ij^i = 2d\jdij — 1 = Xij{Yij^i jYij^2,j ■ ■ ■ ^i-i,j+i)-^t,j+i! (9) 

and the interaction operators in Wen's model take the simple form 

Wi, = (2dyi, - l) (24_,i_^.d,+ij - l) . (10) 

Wen's model transforms differently, depending on the boundary conditions. We 
consider first an N x M lattice with open boundary conditions. Here, the Hamiltonian 
Eq. ([6]) takes the form 

hT"" = ^7 e e (2-^^. - - 1), (11) 

j=i i=i 

where Afij = d\jdij. In this case, neighboring fermions in each row are coupled by 
an Ising-type interactions. However there are no inter-row interactions so that the 
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Hamiltonian takes the form of M — 1 independent Ising chains. With this in mind, it is 
obvious that each row has a two-fold degenerate ground state 

N/2 N/2 

= n 4,1^), 1^"). = n 4-1,1^^), (12) 

i=l i=l 

where \Q) is the state of the vacuum. The states and have anti-ferromagnetic 
excitations in the fermionic bases. The ground state of the Hamiltonian is then given 
by 

M 

m = ^\r^)j, (13) 

i=i 

where n = (rii, n2, ■ ■ ■) and Uj G {+, — }. 

In the case of periodic boundary conditions the diagonalized Hamiltonian takes 
a different form. The periodicity of the system written in terms of spin-half operators 
imposes the constraint Wi+Nj+M = Wij. Then, using the same fermionic transformation 
rules, the Hamiltonian takes the form 

. M N-l M 

^periodic =gY.i: - l)(2^.H-i, - 1) + ^7 E 4> (14) 

j=l i=l j=l 

where Lj is of the form 

N-l N N 

L, = n (2Ar,, - 1) n Ai n (15) 

4=2 i=l i=l 

introducing an interaction between row j — 1 and j + 1 (because the operator affects 
the (z,j — l)-th site). Since the operators Uij and Pij flip the excitation of the fermion 
mode dij, it follows that for even 

N N 

n«.,#^), = m,. n/3^.#^).-i = (le) 

1=1 i=l 

For the sake of simplicity, we deal only with the case that is even in this paper. Then, 
the boundary operators impose that we have to superpose the two eigenstates in 
each row, and one of the ground state of the Hamiltonian f|T^ is given by 

M 

,(i)n 



m') = ®[\r),±\r),). (17) 

The two ground states can be seen as the product of a Greenberger-Horne-Zeilinger type 
entangled state [20] . Further, in the case of periodic boundary conditions and an even 
number of rows (i.e. M even), is not the unique ground state of Wen 's model, 

but we have two other solutions. The complete set of ground state is given by 

M/2 

I^O) = (g) (|V^+)2j ± \r)2,) (|V^+)2,+l ± \r)2,+l) . (18) 

X Note, that if we defined the Hamiltonian Eq. pTjl for periodic boundary conditions, i.e. on a torus, 
the ground state would be 2''^-fold degenerate. However, if we interchanged rows and columns, the 
degeneracy would change to 2^, as we now have a torus of a different size. 
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The ground states of Wen's model can be understood using the stabihzer formahsm 
[2Tj . In the case where g < 0, the ground states of the model will be stabilized by all 
of the terms in the Hamiltonian Wij. In fact, these terms will form a generating set 
for the entire stabilizer group. However, in the case of infinite or periodic boundary 
conditions, the operators Wij do not form an independent set because of the condition 
that riij- W^ij = 1- This implies that Yii^kjj^i^ij = ^ki and if one of the interaction 
operators is removed, then the generating set will become independent. Following the 
stabilizer formalism, if a stabilizer group describing qubits contains k generators in 
the generating set, then the stabilizer group will describe 2^~^ states [21]. As there 
are N operators Wij, of which — 1 are independent, it is clear that in the general 
case periodic lattices will be two-fold degenerate. In the special case that lattice is 
either infinite, or the lattice has dimensions of an even by even number of spins then 
the interaction operators also follow the condition that nj+j=even ^ij = 1- This means 
that there are only N — 2 independent operators in the generating set, so it follows that 
the ground state will be four-fold degenerate. 

5. Connection to Kitaev's Toric Code Model 

In this section, we demonstrate that Wen's model can be transformed into Kitaev's 
toric code model using only local unitary transformations. Since local transformations 
have no influence on topological effects, this allows us to evaluate the entanglement of 
Wen's model by inspecting the ground state of the toric code model. This equivalence, 
coupled with the duality mapping also provides the connection between the two- 
dimensional cluster state and Kitaev's toric code model. The connection between 
Wen's model and Kitaev's toric code model has been studied by Nussinov and Ortiz 
[7j. We investigate this transformation for periodic boundary conditions, and look for 
the necessary conditions for the transformation to be faithful and when the mapping is 
mismatched. Moreover, we show how to construct the exact ground states of Kitaev's 
toric code model in a spin basis. 

The toric code describes a family of simple spin systems with local interactions in 
which the existence of anyons can be demonstrated. Due to their braiding statistics 
the state can been used for topological quantum computation j22]. Kitaev's toric code 
model on a square lattice, 

Hk = -Y.^v-T.Bp, (19) 

V p 

is a sum of constraint operators associated with vertices v and plaquettes p, namely 

A, = l[X„ B, = l[Z„ (20) 

jev jep 

where v = {{i,j), {i + 1, j), (^, j + 1), {i + 1, j + 1)}, Vz, j where i + j = even and 
P = {(^.i). (^ + l.j). (^>J + 1)7 (^ + I7 j + 1)}, Vi,j where i + j = odd. 
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Wen's model can be mapped into the toric code using a completely local unitary 
transformation, UHy^U'^ = where 



Note, that for periodic boundary conditions this relation is only valid for lattices with an 
even number of sites in each direction. Let us briefly explain the action of this unitary: 
If i + j is even, UWtjW = Aij, while for i + j odd we find UWijW = B-ij which is 
a plaquette on p. The conventional method of graphically representing the toric code 
represents the spins as edges on a graph, instead of using the vertices O [6]. One can see 
that this is equivalent by replacing the vertices with diagonal lines. If we then replace 
the spins on even sites with diagonal lines that go from the top left corner to the bottom 
right corner, and spins on odd sites with diagonals that go from the bottom left corner 
to the top right corner then we have a new graph where edges represent spins. This is 
represented in figure [5]on a 4x4 lattice. One can see that, to obtain a faithful mapping 
from Wen's model to the toric code model. Wen's model has to be on a lattice of an 
even number of sites to an even number of sites. If this is not the case then the toric 
code picture will not have alternating diagonal lines across the boundary terms, and 
there will not be perfect star and plaquette operators which consist of only one Pauli 
operator across the boundaries. 



Figure 4. Transformation of a four-by-four square lattice from Wen's spin model into 
Kitaev's toric code model. The blue dots are the actual spins in the four-by-four lattice 
and the yellow dots are the spins which are translated from the first row and the first 
column. The solid line represents the lattice for the Wen model and the dashed lines 
denote the transformed lattice in the toric code model. They are both periodic under 
translation. 




(21) 




\ 
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Having these results at hand, we turn our attention now towards the ground state 
of the toric code. We start with an exphcit construction of the ground states of the toric 
code in the case of periodic boundary conditions. Let Ns denote the number of spins 
on the lattice. One of the four ground states of the toric code can be constructed using 
the set of vertex operators A^, (f = 1, . . . , ^ — 1) as a generating sell§] in the following 
way: 

l^o) = ^^E h ^Tm, (22) 



where n = (ni,n2, . . -ni^ ), Ui G {0, 1} and jfi) denotes the vacuum in the Z-basis. 

2 ^ 

To check that this is indeed the ground state of the toric code, we have to show 
that \ipo) is a stabilizer state of the sets {A^j^, and {Bp}p. We start with the plaquette 
operators. Using the fact that all vertex and all plaquette operators commute, we find 
immediately that 

iXs_i £s_i 

Bp m = E n Ai^Bp \n) = -j:L= E 11 AT \^) = l^o) Vp. 

^2^-1 " k=i V2^~^ " ^==1 

It is slightly more complicated to show that lipo) is also a stabilizer state for the vertex 
operator. To this end, we rewrite the ground state for some fixed v in the following way: 

i^o) = (i+A.)^^E n ^rif^), 

where n„ = (ni, ni, . . . , n^_i, n^+i, . . . , un^ ). Since = 1, we see immediately that 

2 ^ 

A^ lipo) = \ipo) holds. 

The other three ground states I'lpi), ^^2) and \ip3) can be found by applying non- 
contractable loops Wi and W2 composed of Pauli X operators around the torus in the 
horizontal and vertical directions such that the other ground states are {tpi) = Wi IV^o); 
\'^2) = W2 \iIjo) and \ip3) = W1W2 \iJo)- The non-contractible loop operators commute 
with all of the vertex and plaquette operators, thus it is trivial to show that these are 
the other ground states of the toric code Hamiltonian. 



6. Interpreting the Mapping from the Cluster-State to Wen's model on a 
Square Periodic Lattice 

In this section, we investigate the duality mapping between the Cluster state model and 
Wen's model for periodic boundary conditions. Here, special care has to be taken since 
the spins are ordered on the lattice in a diagonal way. For simplicity, we only consider 
the case of an x lattice, where is even. When considering periodic boundary 

§ Note that there arc in fact ^^ star operators because there are two qubits per vertex. However, due 
to the boundary condition J^^^ A^, = 1, only — 1 of them are independent. 
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conditions, we need to be more precise with how we define our indices. We the site 
(A^ + k,j') = {k,j') and the site (1 — k,j') = {N + 1 — k,j'). Then, to ensure the new 
spin operators /i^, fi^ and /i^ obey the SU{2) commutation relations, the mapping has 



is /i^, 



Zi' jv for all i'. Otherwise, all the 



to be changed as follows; when j' 
mappings are the same. 

Having modified the mapping for square periodic lattices, it is then possible 
to map the boundary terms for the cluster state. First, it is easy to show that 
when 2 < j < N — 1, the terms Cij and Cnj map exactly into Wen plaquettes, 

Cij = fi'kj-if^tj-if^Njf^tjy and Cnj = fiN-ij-if^Nj-if^N-ijf^Nj using the standard 
indices. 

The other boundary terms are not as easy to interpret. We have that 

" [i-l 



.k=l 



,n—i+k 
.k=l 



when j = 1 and 



N 



f^i-l,N-lf^tN- 



'N 

.k=0 



X 

i+k,l+k 



.k=l 



N-i+l+k 



when j = N. 

While mapping these boundary terms lead to quite unintuitive terms in the 
Hamiltonian, the ground state of this Hamiltonian can interpreted a little easier by 
manipulating the stablizers. Firstly, one can replace the generators of the stablizer 
group Ci^N with the stablizers Cj_i^iCj^Ar Vi. Then in the Wen picture we have that 
Ci-i^iCi^N = fJ'i-2 i/^f-1 ifJ'i-i N-iK N-iK N- These terms correspond to Wen plaquettes 
interacting between the — 1th row and the 1st row. These plaquettes however are 
skewed, as it is the i — 2th and i — 1th spin on the 1st row that interact with the i — 1th 
and ith spin on the — 1th row. They all also interact with one spin on the A^th row. 
With these terms, we see that we have approximately mapped a cluster state on an 
N X N lattice into a periodic Wen model on an A^ x (A^ — 1) lattice, which interacts 
with the A^th row. It is important to stress however that manipulating the generators of 
the stablizer group only gives us a description of the ground state because if we form a 
Hamiltonian by summing the new generators, there will be a different energy spectrum 
in the excited states. This exchange only allows us to better interpret the ground state. 

Finally, one has to consider the stablizers mapped from the Cj^i terms. In the Wen 
picture, these terms all map into non-contractible loops across the A^ x (A^ — 1) periodic 
lattice in the Wen picture. Then, locally mapping the system into the toric code picture 
reveals that the Cj^i terms map into non-contractible loops of fi^ operators when i is 
odd and non-contractible loops of fi^ operators when i is even. 

It is important that under unitary transformation the energy spectrum is conserved. 
In reference [6], these non-contractible loops are used to store information on a 
topological memory by decreasing the number of available degenerate ground states. 
It is these loops that conserve the degeneracy of the ground state, as we are mapping 
from a cluster-state Hamiltonian which has a non-degenerate ground state. 
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7. Conclusion 

In this work we have reviewed and proven relations between four one-and two- 
dimensional lattice spin models, namely the two-dimensional Cluster State model, Wen's 
model, Kitaev's Toric Code and the one-dimensional Ising model. These mappings can 
be utilized in a plethora of ways. Perhaps the most interesting one is the comparison 
between the two-dimensional cluster state and the Kitaev model. Both models can 
be used as a possible realization of quantum computation by making use of two very 
different ground state properties of the system, i.e. the entanglement and the toplocial 
structure. In this respect, we believe that our work can contribute to new insights 
into questions relating the fields of quantum information science and condensed matter 
physics. For example, can we find a simple way to characterize condensed matter systems 
that can be used for quantum computation? [23] . 

The other important mapping is that between the Wen model and the Ising model. 
Since features of topologically ordered models such as the entanglement structure can 
be difficult to understand intuitively, a mapping between a topologically ordered model 
into something simpler like the Ising model could improve our understanding of the 
behavior of such a system. We believe that an extension of the mapping presented in 
this work to lattices with arbitrary dimension will be a further step to achieve a better 
knowledge of the behavior of topologically ordered systems. 
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